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Abstract 



We propose a mechanism for the observed non-keplerian motion (Tiscareno et al. 2010a) of "propeller" 
moons embedded in Saturn's rings. Our mechanism, in which radial variations in surface density - external 
to, and unaffected by, the embedded moon - result in an equilibrium semimajor axis for the moon due to 
"Type I" angular momentum exchange (Crida et al. 20101, provides a good fit to the observations. Future 
observations should distinguish between our model and others recently proposed. 



1. Introduction 



A massive object embedded in a planetary ring creates a propeller-shaped disturbance in the local 



disk continuum (Spahn and Sremcevic 2000 Sremcevic et al. 2002 Seifi et al. 2005). Swarms of such 



"propellers," each with central moonlet radii r ~ 100 m, occur in the 3,000-km-wide "Propeller Belts" 



region in the middle part of Saturn's A ring (Tiscareno et al. 2006 2008 Sremcevic et al. 20071, though 



photometric ambiguities make it difficult to precisely ascertain the central moonlets' masses (Tiscareno 



et al|[2010bj. 

A second propeller-rich region has been identified in the outermost regions of the A ring ( |Tiscareno et al. 



2010a| hereafter T10). Propellers in this region are both larger and rarer than in the Propeller Belts, and 



several have been tracked for periods of several years. The propeller nicknamed "Bleriot'j^] is the largest 
(r ~ 1 km) and best-tracked (> 100 detections over 5 years). Its orbit is predominantly keplerian, with 
longitude residuals less than ±0.15° (200 km) over a period of nearly 5 yr. However, the residuals are much 
larger than measurement error and are clearly systematic (T10; data reproduced in Fig. [I]). If the observed 
non-keplerian motion is interpreted in terms of changes in the instantaneous semimajor axis a, then Bleriot 
migrated outward from mid-2006 to mid-2007 at a rate of a — +0.11 km/yr, and inward from late-2007 to 
early-2009 at a rate of a = -0.04 km/yr (T10). 

A linear trend in semimajor axis, that is to say constant a, corresponds to a quadratic trend in the 
longitude residual, X(t) oc t 2 . The observed longitude residuals for Bleriot can be interpreted in terms of a 
series of piecewise quadratics (red lines in Fig. [T]) that would imply episodic migration, perhaps caused and 
punctuated by periodic encounters or collisions ( |Lewis and Stewart 2009 Kirsh et al. 2009, T10). While 
plausible, this mechanism needs further development, and this paper will not discuss it further. Other 
mechanisms for Bleriot 's non-keplerian motion have been proposed, and are reviewed below: 

The observed longitude residuals for Bleriot can be interpreted in terms of a sinusoidal oscillation with 
period 3.7 yr (blue line in Fig. [TJ. Although the most obvious physical mechanism to produce a sinusoidal 
longitude residual is resonant interaction with a larger moon exterior to the rings, a mechanism which 
governs many phenomena in the rings (e.g., |Tiscareno| 2012), no resonance is known near Bleriot's position 
that could plausibly cause longitude residuals with the observed large amplitudes. In the "frog" mechanism 



suggested by Pan and Chiang (2010 hereafter PC10) and further developed by Pan and Chiang ( 2012[ ), the 
propeller moonlet interacts primarily with the mass at either end of the propeller gap. Modeling those ends 



1 Propellers that have been observed repeatedly are nicknamed to facilitate identification. 
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Figure 1: Observed longitude of the propeller "Bleriot" over 4 years, with linear trend Ao = 616.7819329° day -1 subtracted 
off. Data reproduced from T10, q.v. for more detail. Panel (a) contains all the data, while panels (b), (c), and (d) contain 
subsets of the data shown in greater detail. The blue line indicates a linear-plus-sinusoidal fit to all the data, while the red lines 
indicate piecewise quadratic fits corresponding to a constant drift in semimajor axis (T10). The black lines indicate exponential 
fits as described in Section [2] The data from 2005 and 2009 are too sparse to distinguish among fitting functions, so are fit to 
a simple linear trend. 



as co-orbiting masses, PC10 found that the propeller moonlet plausibly librates with the observed amplitude 
and period. However, one wonders whether the moon-formed gap responds sluggishly enough to allow the 
moon to librate within it. Two predictions of this may soon be tested: that Bleriot's 3,000-km-long gap 
structure (which is only seen clearly in a few images) is stationary with respect to the librations, and that 
the longitude residual will continue to follow a sinusoidal profile into the future. 

Other hypotheses for Bleriot's non-keplerian motion, including the one we present below, rely on the 
concept of "Type I migration." As classically formulated for protoplanetary disks (Ward 198G 1997 



Papaloizou et al. 2007), the angular momentum exchange at inner Lindblad resonances between a disk and 



an embedded mass fails to exactly cancel with that at outer Lindblad resonances, resulting in a differential 
torque that leads to inward migration of the embedded mass. However, classical Type I migration depends 
crucially on the gas component of the disk, which causes Lindblad resonance locations to shift asymmetrically. 
For the case of planetary rings, which are strictly particulate, Crida et al. (2010, hereafter C10) re-derived the 
equations for Type I migration from first principles, using analytical arguments and numerical simulations to 
trace the angular momentum exchange between streamlines of continuum ring particles and the embedded 
moon, and deriving a profile of torque per unit disk surface density as a function of the impact parameter 
b (i.e., the difference in semimajor axis) of ring particle streamlines with respect to the moon (Fig. [2]). The 
primary contribution does not depend on resonances, which are more symmetrically placed in a non-gaseous 
particulate disk, but on intrinsic asymmetries in the single-impulse transfer of angular momentum. 

C10 find an intrinsic asymmetric torque, akin to classic Type I migration, in the case of a particulate 
disk with uniform surface density; it is always inward (while Bleriot is seen to move both inward and 
outward) and is one to two orders of magnitude too weak to explain the magnitude of Bleriot's observed 
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Figure 2: Angular momentum exchange AJ between a ring streamline and an embedded moon as a function of the normalized 
impact parameter b = b/rjj. AJ is normalized by the moon's specific angular momentum J m = a^no and by its Hill radius 
r H . Data from |Crida et alT] | |2010| l. 



non-keplerian motion (CIO). However, it is known that Saturn's rings do not have uniform surface density. 
Rein and Papaloizou ( 2010| hereafter RP10) considered stochastic temporal variations of the surface density 



due to self-gravity wakes in the region of maximum angular momentum transfer, finding that they induce 
a random walk in the propeller moon's semimajor axis. Under certain realistic conditions, the magnitude 
of the longitude residual approaches that observed for Bleriot, and since a random walk has no preferred 
frequency, any frequency might be apparent over a limited period of time. The model of RP10 predicts 
that the quasi-sinusoidal behavior observed by T10 will not continue, but will instead be seen as part of an 
overall random progression. 

2. Analytical Model 

In Saturn's rings, the surface density changes not only with time (as in the model of RP10), but also 
with position (ring radius, or the distance from Saturn's center within the ring plane). Thus, we propose 
here that Type I torques can act to keep the embedded moon's semimajor axis at an equilibrium position, 
given a certain radial profile of the surface density. Indeed, radial surface density structure in the ring means 
that any propeller must constantly migrate due to imbalanced torques, unless it finds a location where the 
torques balance. Thus, we should expect to find them at such equilibrium points. 

The "sweet spots"[^] where angular momentum exchange |AJ(6)| is maximized in the CIO calculations of 
Type I torque are at ~ 2 Hill radij^jon either side of the embedded moon's semimajor axis (Fig. [2]). These 
local maxima occur because AJ(b) increases with impact parameter b in the horseshoe region (b < 1.8rjj), 
while it decreases with b in the farther region (b > 2.5r#) where particles are only moderately deflected by 
the encounter (see, e.g., Fig. 3.30 of |Murray and Dermott| [T999 ). Between those two regions lies a chaotic 



zone where AJ(b) is depressed by the fact that particles can be deflected to either side of the moon, thus 
leading to two peaks at b ~ 1.8r# and b ~ 2.5r# (CIO). The sign of the torque is positive (negative) on the 
inward (outward) side of the moon in the moderate deflection region, and vice-versa for the horseshoe and 
chaotic regions (Fig. [2]). 



In baseball, the "sweet spot" is the location on the bat at which a struck ball is most efficiently propelled. So, here, we refer to 
the impact parameter at which angular momentum is most efficiently transferred between the moon and the disk. 

^An object's Hill radius th — ao(m/3Af) 1//3 , using variables defined later in this section, is the approximate extent of its gravita- 
tional dominance. 
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Although ClO's results are plotted for only a single value of the embedded moon mass, scaling relations 
can be determined. Such scaling is needed because the moon mass used in their models is 400 times larger 



than Bleriot's likely mass (T10) - in fact, it is 20 times larger than the moon Daphnis' mass (Porco et al 



2007) and thus represents a moon that would not create a propeller structure but rather a fully cleared gap. 
C10 state that their total cumulative torque scales as m 4 / 3 oc r^, where m is the mass of the embedded moon. 
However]^] one factor of m 1 / 3 oc rn does not enter into our calculations since it comes from integrating the 
difference between the torques on the inner and outer sides of the moon under ClO's assumption of constant 
surface density, as their SJ(b) = AJ(b) — AJ(—b) oc AJ(b) • bruja-a (Eq. 34 of C10), using the scaled impact 
parameter b = b/rn- For our purposes, we note that AJ(b) itself is proportionaj^jto J m r#, where the moon's 
specific angular momentum is J m = a^no for semimajor axis ao and mean motion no (for this reason, in 
Fig. [ij we normalize AJ(b) by rn as well as by J m , obtaining a curve that applies to embedded moons of 
any mass), and that a further factor of r 2 H arises when calculating the integrated torque (rewriting Eq. 36 
of C10): 

T = 3n r 2 H / a(b')A J(b')b'db'. (1) 



Thus, for our case in which a(b) must be convolved with the angular momentum exchange AJ(b), rather 
than being held as a constant so that SJ(b) becomes important, the total torque T oc r^ H oc m. Since 
migration rate a oc T/m, we conclude that a given variation in surface density will induce the same a in 
moons of different masses, though the radial scale over which the surface density variation occurs must scale 
with th - so that a(b) remains constant - to preserve the effect. 

Consider a ring whose surface density varies radially but is azimuthally symmetric, and which has a local 
minimum in surface density at a given radius, as in Fig. [3] The radial density structure a(b) is external to, 
and unaffected by, the embedded moon; in fact, such radial density structure is common in Saturn's rings 
(see Section [3] for more discussion). The moon's evolution is primarily driven by variations in a(b) at the 
"sweet spots" on either side of the moon. The total torque T on the moon is found by integrating Eq. [l] 
convolving AJ(6) with cr(6); a stable equilibrium can be created simply by setting dcr/d6 negative on the 
inward side of the moon and positive on the outward side. If the embedded moon has its semimajor axis 
suddenly altered from this equilibrium position by an encounter (ILewis and Stewart I 2009), for example 



with a large fluctuation of self-gravity wakes or with a large ring particle, the torque from the surrounding 
material will then be out of balance and will push the moon's semimajor axis back toward the equilibrium 
position (Fig. |3|. 

In the following paragraphs we analytically consider the simplest case, in which da/db has a constant 
negative value at the inward sweet spot, and the same constant value (except positive) at the outward sweet 
spot, and the equilibrium semimajor axis ao lies exactly between the two regions. If we simplify the angular 
momentum exchange profile as two delta functions, one on each side of the moon, 

— = AS(b±b ), (2) 

G 

where A and 60 are free parameters (we use A = 4.7 x 10 18 cm 4 s~ 2 , a value obtained by integrating Eq. [l] 
with AJ(b) given by Fig. [2] and aib) set equal to unity; and 60 = 2.78 km in order to mimic the balance 
point of the full exchange function), then the convolved torque is given by 

T = -2A^-(a-a ) 1 (3) 



The insight of A. Crida (personal communication, 2011) was invaluable in formulating the analysis given in the remainder of this 
paragraph. 

^The asymptotic relations shown in Fig. 3 of C10 make this clear. These expressions can be written as A J(b) ~ (brn /a>o)Jm for the 
horseshoe region and A.J(b) ™ {brn /2ao) Jm for the chaotic deflection region. For the moderate deflection region, AJ(b) as expressed 
in Eq. 32 of C10 is proportional to m 2 /b 5 oc r/f/6 5 . The inference that AJ(6) oc rn is confirmed in direct calculations by A. Crida 
(personal communication, 2011). 
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Table 1: Fit parameters to A — n t = p 1 + (po - p 1 ) e - 2 ' r ( t -P3)/p2 f or data in Fig. FT 



Year 


Po 


Pi 


P2 


Pa 


mid-2006 to carly-2007 


0.00 


0.12 


0.98 


2006.59 


lato-2007 to carly-2009 


0.066 


-0.098 


3.21 


2007.94 



where ao is the equilibrium semimajor axis and a is the moon's actual semimajor axis. The angular momen- 
tum is given by L = m(GMa) 1 ^ 2 , assuming that any non-zero eccentricity will be quickly damped away, 
where G is Newton's constant and M and m are the masses of Saturn and the moonlet. We take the latter 
as m ~ 2 x 10 15 g ~ 4 x 10~ 15 M for an icy moon with radius ~ 1 km (T10) and internal density close to 



the Roche critical density p c « 0.5 g cm d (Porco et al. 2007). The torque is then T = L = mnoaoa/2 



Setting this equal to Eq. [3j we obtain a simple differential equation for the semimajor axis offset a — a : 

d(a - ao) 



dt 

where 



+ cj (a-a ) =0, (4) 



4^4 da 

^o = (5) 

mnofflo dr 

The solution, with initial displacement a = ad at t = 0, is 

(a-a ) = (a d -a )e- u ' ot . (6) 

Using the standard equation for keplerian shear and substituting Eq. [6j the longitude residual A — Ao, where 
Ao = n Q t, will accumulate as 

^r^ = -l^(a d -a )e-^. (7) 
at z ao 

Integrating, we obtain 

\-\ = l^(a d -a )(e-»° t -l). (8) 
2 a uj 

Note that, as t — > oo, the longitude residual asymptotically approaches a constant negative (positive) value 
for a positive (negative) initial offset in semimajor axis. 

This simplified model result can be compared to data. In Fig. [I] the observed longitude residuals in two 
well-sampled periods are separately fit to functions with the form of Eq. [8] (see black lines in Fig. [I] and 
Table[I]). In the interval from mid-2006 to early-2007 (Fig.flJ;), we obtained (2tt/luq) = 1-0 yr with A — Ao = 
270 km = 0.0020 radians. Thus, taking n = 616.782° day" 1 at a = 134912 km (T10), we calculate from 
Eq.[5]that da/db = 3.8 g cm~ 2 km -1 and from Eq.[8]that — ao = —300 m. Similarly, in the interval from 
late-2007 to early-2009 (Fig. [TJi), we obtained {2tt/uj ) = 3.2 yr with A - A = -390 km = -0.0029 radians, 
yielding der/d6 = 1.2 g cm -2 km -1 and — ao = 130 m. 

The initiation of our mechanism is not simply explained. Our model requires three events (where the 
black lines cross in Fig. [lj over our five-year data set at which there is a rapid change in a on the order 
of a few x 100 m. Direct delivery of sufficient momentum could occur if an object with approximately one- 
thirtieth the mass of Bleriot (i.e., one-third its radius) strikes at Bleriot's escape velocity of ~ 50 cm s _1 , 
The size distribution reported by T10 predicts that approximately one such object should occur per annulus 
of radial width 35 km, while the radial cross-section for interacting with Bleriot is at least the diameter of 
its central moonlet ~ 2 km and more likely the total radial width of the cleared region ~ 15 km - this is the 
sum of the radial offset between the two lobes Ar ~ 5 km (T10) and the radial full-width at half-maximum 
of a lobe, which we estimate from the T10 dataset to be ~ 10 km. Another possibility is interaction with an 



unusually dense (~once-per-year) self-gravity wake (SGW; see Cuzzi et al. (2010)) of similar mass. Little is 
known about the detailed dynamics of SGWs, especially infrequent outliers. Any resulting eccentricity need 
not be quickly damped; there is no empirical evidence to discount a small eccentricity ~ 10~ 6 for Bleriot, 
which would lead to azimuthal residuals ~ 0.1 km. 
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Figure 3: In this toy model, the embedded moon's semimajor axis (dotted blue line) begins with an initial inward displacement 
from equilibrium at t = (a) and trends back to its equilibrium value (b). The blue and cyan curves are the angular momentum 
exchange profile (positive and negative contributions, respectively), here approximated as a pair of narrow gaussians (Eq.[2]l. 
The black curve is an input radial surface density (a) profile, with dcr/db = —4.8 g cm -2 km -1 on the inward side of the 
equilibrium semimajor axis and +1.44 g cm -2 km -1 on the outward side. Both slopes occur only in the region of maximum 
angular momentum exchange. The right-hand panels show (c) the semimajor axis offset Aa, (d) the azimuthal residual AA, 
and (e) the torque T as functions of time. 



We also note that the method of Rein and Papaloizou (2010), which has kicks from SGWs at its core 
while we mention such kicks above only as a possible initiation for our mechanism that focuses on disk-moon 
torques, similarly suffers from problems with the magnitude of the available perturbations. 

However, we believe this model remains worthy of consideration because the data (Fig. [T]) resemble a 
series of curves somewhat better than they resemble a smooth sinusoid. The residuals for this new functional 
fit are significantly better than for the sinusoidal fit (though with 10 degrees of freedom, compared with 5 
for the sinusoidal fit), and are comparable to those for the piecewise quadratic fit (T10), with particular 
improvement for the spring 2007 data (Fig. [TJ:) . 

The present model requires only three "kick" events to explain the 5-year data set, since the 2006/07 and 
the 2007/09 lines intersect each other. However, the differing values of ojq for the two segments may indicate 
that da/db does not have the same value at the inward and outward "sweet spots", leading to differing 
response times to an offset. With this in mind, we construct a simple numerical version of our model in the 
next section, using the entire angular momentum exchange profile rather than reducing it to a pair of delta 
functions as we did for the above analytical treatment. 



3. Numerical Model 

We verified the above analytical results by numerically integrating the C10 equations. 

Firstly, we represented the angular momentum exchange function as a pair of narrow gaussians after 
the manner of Eq. [2] In two separate runs, with a radial surface density profile consisting of a "valley" 
with slopes of dcr/d& = 3.8 (1.2) g cm~ 2 km -1 in the vicinity of bo, and with an initial semimajor axis 
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Figure 4: The same as Fig. [3] but with the full angular momentum exchange profile (the profile shown in Fig. [2] multipled by 



the surface density profile), and with a double-troughed surface-density profile with da/db = ± 13.5 g cm 2 km 
of maximum angular momentum exchange. 



in regions 



—300 (+130) 111 from the e, we indeed found that the semimiajor axis took about 1.0 (3.2) yr to return to 
equilibrium and that the total azimuthal residual came to about +270 (—390) km. 

Secondly, to produce a model that does not invoke temporal changes in the surface density profile, we 
used a profile with a steeper slope on the inward side and a gentler slope on the outward side (Figs. [3] and [5]). 
In this way, only the inward (outward) slope is encountered by the exchange functions when the semimajor 
axis is offset inward (outward) , and different response times for the the two cases are thus possible. We again 
obtained results similar to the observations when the "valley" floor was 3.8 g cm -2 below the background. 

Finally, we extend our treatment to the case of the actual angular momentum exchange function (Fig. [2]). 
There are many ways to configure the surface density so that the longitude residual matches the data, 
including a square well and the modified square well seen in Figs. [| and |5j but we found the most efficient 
model (requiring the smallest Act) to be that shown in Figs. [4] and [6j Because the angular momentum 
exchange is spread over a larger annulus, it is harder to have steeper slopes encountered only for an inward 
initial offset and gentler slopes only for an outward initial offset, as in the previous paragraph. In the 
example shown, the "valleys" are 5.4 (1.42) g cm -2 lower than the background surface density for initial 
semimajor axis offset in the inward (outward) direction, yielding slopes of da/db = 13.5 (3.55) g cm~ 2 km~ 2 . 



4. Discussion and Conclusions 

We have shown that a modest permanent radial variation in the surface density, combined with plausible 
occasional "kicks" to the semimajor axis, can lead to non-keplerian motions of an embedded moonlet that 
are similar to those observed for Bleriot. It finally remains to inquire whether such variations may actually 
exist in Bleriot's vicinity. This region is characterized by variable structure (Fig. [7]), most likely due to spiral 
density waves. Our model is viable if, when averaged over orbital timescales (which are short compared 
to the migration rate), the mean torques due to that structure yield anomalies on either side of Bleriot's 
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Figure 5: This figure is identical to Fig. [3] but for an outward initial displacement. The vertical dashed line in panels (c), (d), 
and (e) indicates the approximate length of the observation interval (Fig. lip). 



position. Optical-depth variations in this part of the A ring are ±0.2 (M. Hedman, personal communication, 
2011). Using a mass extinction coefficient for the outer-A ring of 0.01-0.02 cm 2 g _1 (Colwell et al. 2009), 
this corresponds to surface density variations as high as ±10 g cm" 2 . The mechanism described herein 
requires surface density anomalies of only a few g cm~ 2 . 

Although we cannot prove whether the observed surface density variations are underlain by the much 
smaller average anomalies necessary to activate our mechanism, nor is such proof likely to be possible 




8 




0.96 L , , , , , , , . . . , , . ; . , , . . . . , , , , , , . 

134900 134905 134910 134915 134920 134925 

Radius (km) 



Figure 7: Normalized brightness (I/F) profiles of the A ring in the vicinity of Bleriot's semimajor axis (vertical dotted line) 
from the following images: N1540681631, taken on 2006 October 27 (green); N1541717040, taken on 2006 November 8 (cyan); 
N1560309735, taken on 2007 June 12 (blue); N1595337202, taken on 2008 July 21 (red); and N1654254507, taken on 2010 June 3 
(black). Regions of maximum angular momentum exchange are shaded in gray. 



without greatly improved observations and/or difficult numerical modeling, our purpose here is only to 
argue that it is plausible. The model presented in this paper is motivated by the observed orbital evolution 
of Bleriot and other propellers, which seems to be characterized by periodic "kicks," with slow variations 
in the semimajor axis (seen as curvature in the longitude residuals) occurring between those kicks. It is 
in response to these observations that we have constructed a mechanism that is dynamically plausible and 
consistent with existing data. 

Our proposed mechanism joins two others recently outlined, and all three mechanisms make contrast- 
ing predictions regarding the future behavior of Bleriot's azimuthal residual. While PC10 predict that the 
quasi-sinusoidal trend discerned in the azimuthal residual to date will continue, and RP10 predict that 
variations in the azimuthal residual will be seen in the future to result from stochastic variation in the 
semimajor axis, we predict periodic (less than one per year, on average, from our interpretation of existing 
data) sudden "kicks" to the semimajor axis followed by asymptotic return to the equilibrium value. Future 
data, beginning in 2013 when Cassini will again spend significant time out of Saturn's equatorial plane with 
good viewing of the rings, should distinguish among these models. 
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